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Abstract 



Rare B meson decays are an excellent probe for beyond the Standard 
Model physics. Two very sensitive processes are the B® d — > and 
b — > sfi + fi" decays. We report recent results at a center of mass energy of 
a/s = 1.96 TeV from the CDF II and DO collaborations using between 3.7 
fb" 1 and 6.9 fb _1 taken during Run II of the Fermilab Tevatron Collider. 
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1 b — > sfi + fi decays 

b — > sfi + fi~ decays are flavor changing neutral current (FCNC) processes that can 
only occur through higher order box or penguin amplitudes in the Standard Model. 
New physics can be probed by measuring various combinations of their decay rates. 
One of the most sensitive observables is the forward-backward asymmetry of the 
muons (Afb) as a function of the squared di-muon mass. 

Three decays of interest are B + B° d //+ pT A"*(892)° , and 5° -> 

fi + fi~4>. The first two decays were observed at BaBar, Belle, and CDF. Belle found a 
2.7a deviation for Aps for B® — > l + l~K* pp. We describe hre the CDF analysis that 
uses 4.4 fb _1 of data 0. 

1.1 Analysis Method 

1.1.1 Branching Ratio Measurement 

Branching ratios for B G d -»■ /j, + fi~K* (892)°, B + ->■ /j, + fi~K + , and B® -> are 
measured relative to normalization modes, where the two muons originate from a 
J/V decay. For the event reconstruction CDF requires two muons with a transverse 
momentum (p?) greater than either 1.5 GeV/c or 2.0 GeV/c depending on the muon 
trigger. The three modes are then reconstructed where the K*(892)° is reconstructed 
from K*(892)° — > K + n~ and the <j) is reconstructed from <fi — > K + K~. To avoid 
contamination from resonant modes such as the J/ip and i/)', candidates with di-muon 
masses near these resonances are rejected. 

The events then have to meet loose preselection requirements before an artificial 
neural network (NN), which combines multiple discriminating variables into one vari- 
able, is applied. Signal is modeled with py-reweighted Pythia signal Monte Carlo 
simulations (MC). The reweighing is done by comparing the MC pr distribution with 
that of the normalization modes. The background is modeled by sampling the B me- 
son mass sideband regions. For the B® — > fi + fi~K* (892)° and B + — > fi + fi~K + modes 
there are significant physics backgrounds in the lower sideband region and thus only 
the higher sideband region between 5a and 15a (er=20MeV/c 2 ) above the known B 
meson masses is used. For B® — > p, + pT4> both lower and upper sidebands are used. 

The final signal yield is obtained by an unbinned maximum likelihood fit to the B 
meson mass distribution. The probability distribution function (PDF) of the signal is 
parametrized with two Gaussian with different means while the the background PDF 
is described by a first or second order polynomial. Peaking background contributions 
are subtracted from the fit results for the signal yields. The only significant peaking 
contribution is cross-talk among B d — > n + n~K* (892)° and B® — > which has a 

~1% contribution to the total observed signal MC yields. The final branching ratios 
are calculated as follows: 
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jyNN gloose 

B(B -> fc„V) = -i^T ■ J£^-M^- ■ B(Jm ■ B(J/i> -+ M V), (1) 
where B signifies the B% B + , or B° s and /i represents if* (892)°, A+, or 0, 

^loosc 

and Njl^ft are the yields after the optimal NN, e iooit h is the relative efficiency of the 

loose selection cuts, and eF£+ _ is the NN cut efficiency on the loosely-selected events. 
The NN is not applied to the normalization mode because the signal/purity and size 
is sufficient with the loose selection cuts. The NN cut efficiency are obtained from 
signal MC. 

The three leading systematics are the systematics on the efficiency, B(J/iph), and 
background PDF. The main sources within the efficiency systematic errors are the 
MC reweighing and the NN cut. The total systematics for the B® ->■ A'* (892)°, 
B + — > fi + fi~K + , and B® — > are 9%, 7%, and 32%, respectively. 



1.1.2 Forward-Backward Asymmetry Measurement 

For the B® — > fi + fi~ K* (892)° decay the forward-backward asymmetry (Aps) as well 
as the K* longitudinal polarization (Fl) are measured. These are extracted from 
cos# M , the cosine of the helicity angle between the /i + momentum vector and 
the opposite of B (B) meson momentum vector in the di-muon rest frame, and cosOk, 
the cosine of the angle between kaon momentum and the opposite of the B meson 
momentum vector in the K* (892)° rest frame. B decay amplitudes are calculated 
using operator product expansion and Wilson coefficients. There are many non-SM 
predictions for Apb from models with different Wilson coefficients [3]. 

Afb and Fl are extracted using a unbinned maximum likelihood fit containing B 
mass shape signal and background PDF's as well as signal and background angular 
shape PDF's. The mass shape PDF's are divided into several di-muon mass bins 
and are described as 25 bin histograms. The combinatorial background PDF is taken 
from the B meson higher sideband. The angular acceptances are also described as 25 
bin histograms and are derived from phase space signal MC. 

As a control A FB and F L are fitted to B% -t J/^A*(892)° and A FB only to 
B + — > J/ipK + . This cross check yielded measurements that were consistent with 
other measurements. 

The systematic uncertainties for the Fl measurement in B^ — > fi + fi~ A*(892)° 
are dominated by a fit bias ranging between 2% and 7% for the different di-muon 
mass bins. This bias affects teh estimates of Afb or Fl when the true value of 
these quantities are close to physical boundaries. The lowest di-muon mass bin also 
has a significant contribution from the angular efficiency systematic (~5%) attained 
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from changing the binning of the angular acceptance histograms. The total of all the 
systematics ranges between 2% and 8% for the different di-muon mass bins. 

For the Afb measurement in B d — > fi + fi~ K*(892)° the dominant systematic un- 
certainties are the angular efficiency, fit bias, and the statistical uncertainty from the 
Fl fit. The Fl fit systematic comes from varying Fl by ±ler where a is the statistical 
error from the Fl fit. The total systematic errors range is 1%-10%. 

The most significant systematic uncertainties in the Aps measurement for B + — > 
fi + fi~K + is given by the angular efficiency and angular background. Angular back- 
ground is estimated from the B higher mass sideband. The systematic is estimated 
by changing the sideband region. The total systematic uncertainty ranges between 
3% and 8% for the different di-muon mass bins. 



1.2 Results 

The resulting yield for the three decays are shown in Figure [U All measured branching 
ratios agree with previous measured values as well as theoretical predictions. CDF re- 
ports the first observation of B® — > fi + fi~ 0with a significance of ~ 6a and a measured 
branching ratio of B(B° S fi + fi~<j)) = (1.44 ± 0.33[stat]±0.46[syst]) x 1(T 6 . 




Figure 1: Yields for B° d -> fi+ \r K* (892)° , B + ->■ fi + fi~K+, and B° s ->■ from 
left to right. In all histograms the red vertical lines indicate the signal regions. 

The results of the polarization and forward-backward asymmetry measurement 
shown in Figure [2] (Aps f° r B + — > fi + fi~K + is shown as a cross check) may contribute 
to contrain several beyond SM physics models. They are compatible and competitive 
with the results from the B- factories [I], jl], [5]. 



3 




q 2 (GeV 2 /c 2 ) q 2 (GeV 2 /c 2 ) q 2 (GeV 2 /c 2 ) 



Figure 2: K* Longitudinal Polarization, A FB for B° d — > fi + fi~ K* (892)°, and A FB 
for B + — > fi + fi~K + (from left to right). The green hatched areas are due to the 
charmonium veto (J/ip, ift')- 

2 B® d — > fi + fi~ decays 

B° s d — > fi + fi~ are FCNC decays that are highly suppressed by the SM. The SM 
predictions for these branching fractions are B(B® — > = (3.2 ± 0.2) x 10~ 9 

and B(B d — > /i + /i~) = (1.00 ±0.1) x 10~ 10 [5J. These predictions are one order 
of magnitude smaller than the current experimental sensitivity. Previous bounds 
from the CDF collaboration, based on 3.7 fb" 1 of integrated luminosity, are B(B® — > 
/i+/i") < 4.3 x 10" 8 and B(B° d -» fi+ >-) < 7.6 x 10" 9 at 95% C.L. 0. DO bounds were 
previously set at B(B° S < 5.1 x 10" 8 at 95% C.L. using 6.1 fb" 1 of data [8]. 

LHCb bounds of B(B° S < 5.6 x 10~ 8 and B(B° S /i+/i") < 15 x 10" 9 at 

95% C.L were recently set using 37 pb -1 of data 

Enhancements to the rate of B® — > fi + fi~ decays occur in a variety of different 
new-physics models. In supersymmetric (SUSY) models, new particles can increase 
B(B® — > fi + fi~) by several orders of magnitude at large tan/3. In the minimal su- 
persymmetric standard model (MSSM), the enhancement is proportional to tan 6 /3. 
Global analysis including all existing experimental constraints suggest that the large 
tan/3 region is of interest [IUJ [TTJ [T2]. For large tan /3, this search is one of the most 
sensitive probes of new physics available at the Tevatron experiments. In contrast, 
SUSY R-parity violating models [10] and non-minimal flavor violating models [13] 
can both enhance B® — > and B d — > separately even at low tan/3. In the 

absence of an observation, limits on B(B° S — > /i + /i~) are complementary to those pro- 
vided by other experimental measurements, and together would constrain the allowed 
supersymmetric parameter space. 

This document describes current status from DO and CDF, as well as the planned 
updated analysis by CDF. 
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2.1 Analysis Method 



Both CDF and DO collect opposite sign muon candidate events using di-muon triggers. 

B + — > J/ipK + events are collected on the same triggers as a relative normalization 
mode to estimate £>(-£>!? — > /i + yU~) as: 

B(B° -+ „ V) = • • £ • B{B+ -+ J/^K+) ■ B{J/i> -+ „ V), (2) 

where Nbo is the number of candidate B® — > events, a^o is the geometric 

and kinematic acceptance of the di-muon trigger for B® — > ji + ji~ decays, e^o al is the 
total efficiency (including trigger, reconstruction and analysis requirements) for B® — > 
events in the acceptance, with Ng+, and e^+ al similarly defined for B + — > 
J/i/jK + decays. The ratio f u /f s accounts for the different 6-quark fragmentation 
probabilities and is (0.402 ± 0.013)/(0.112 ± 0.013) = 3.589 ± 0.374, including the 
(anti-)correlation between the uncertainties [H]. The final two terms are the relevant 
branching ratios B(B+ -> J/ifjK+) ■ B(J/i/j ->■ fi + fi~) = (1.01 ±0.03) x 10~ 3 • (5.93 ± 
0.06) x 10~ 2 = (6.01 ± 0.21) x 10~ 5 [II]. 

The CDF analysis described is also sensitive to — > fi + fi~ decays. The value of 
B(B® — > /i + /i~) is estimated from Equation [2] substituting B® for B®, and changing 
fu/fs f° ful fd — 1- All other aspects are the same as the B° s — > /i + fi~ search 
except where noted below. The analysis is done by first estimating the acceptances 

and efficiencies, then creating a multivariate discriminant for background rejection. 
This discriminant is optimized with Pythia signal MC and data mass sideband and 
validated with the B + — > J/ipK + normalization. The background is then estimated, 
which has two sources: combinatorial background and peaking background (B— > 
h + h~). Finally when the background is well understood the signal region is unblinded. 



2.2 Signal and Background Properties 

The signal candidates are fully reconstructed events with a secondary vertex due to 
the long lifetime (~ 450/im) of the B® meson. Signal events will have a primary- 
to-secondary vertex vector that is aligned with the B® candidate momentum vector 
(Figure |3]). Another property of signal events that is unique is that they are very 
isolated (with few tracks near the muon tracks) due to the hard B fragmentation. 

Background events tend to be partially reconstructed and be shorter lived than 
signal. They also have a softer pt spectrum, higher activity of tracks, and misaligned 
primary-to-secondary vertex and momentum vectors. The combinatorial background 
consists of sequential semi-leptonic decay (b — > cfi~X — > fi + fi~X) and and double 
semi-leptonic decay {bb — > /i~/i + X). The expected number of combinatorial back- 
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Figure 3: Sketch of the decay topology. 



ground events in the signal window is estimated by extrapolating the number of 
events in the sideband regions to the signal window using a first order polynomial fit. 

The peaking background from two-body hadronic B decays is also evaluated. 
These backgrounds are about a factor ten smaller than the combinatorial background 
and need to be estimated separately. The dominant sources of background are the 
decays of B° s and B° to h + h'~ final states, where h or h' can be either n^ 1 or and 
are misreconstructed (fake muons). 

2.3 CDF Analysis with 3.7 fb 1 of Data 
2.3.1 Signal Discrimination 

Six variables are used for the signal discrimination (Figure H]). Three of the variables 
are secondary vertex variables. These variables are combined in an artificial neural 
network which is optimized with signal MC and data mass sidebands (Figure HI left). 
The NN output is then divided into three bins while the di-muon mass is split into 5 
mass bins for limit calculation. 

A: the proper decay length; 

Aa: the 3D opening angle between the B® flight direction, p(B), and the direction 
of the decay vertex - estimated as the vector originating at the primary vertex 
and terminating at the muon-pair vertex; 

Iso: the isolation of the candidate B" defined as Iso = p T a /{p T s + J2iPt(AR < 1.0)), 
where the sum is over all tracks within an rj — <fi cone of radius R = 1.0, centered 
on p(B); 

\/<j\: proper decay length significance (A divided by the uncertainty on A); 

Pt(B°): the transverse momentum of the B° s candidate; 

Pt(/x): the transverse momentum of lower momentum daughter muon. 
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Figure 4: Left: Neural network output for signal MC and data mass sideband. 
Right: Distributions of neural network input variable for signal MC and data mass 
sideband. 

2.3.2 Control Samples for Combinatorial Background Estimates 

CDF cross checked the background estimate procedure on four independent back- 
ground samples: 

OS-: opposite-sign muon pairs, which pass the baseline and vertex cuts and have a 
negative lifetime; 

SS+: same-sign muon pairs, which pass looser baseline and vertex cuts and have a 
positive lifetime; 

SS-: same-sign muon pairs, which pass looser baseline and vertex cuts and have a 
negative lifetime; 

FM+- : opposite-sign fake-muon pairs, at least one leg of which is required to fail the 
muon ID requirement, passing looser baseline and vertex cuts for both positive 
and negative lifetimes. 

These are representative of various background contributions and thus are ex- 
cellent control regions to test the method of background estimation. We compare 
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the estimated background with actual observation. The result of this cross check is 
summarized in Table [TJ The agreement between predicted and observed number of 
background events in the signal region is good across all of our control regions. The 
backgrounds were estimated using the same method as for the signal sample. 







CMU-CMU 


CMU-CMX 


sample 


NN cut 


pred 


obsv 


prob(%) 


pred 


obsv 


prob(%) 




80 < vmm < 95 


275 ± (9) 


287 


26 


310 ± (10) 


304 


39 


OS- 


0.95 < v NN < 0.995 


122 ± (6) 


121 


46 


124 ± (6) 


148 


3.2 




0.995 < v NN < 1.0 


44 ± (4) 


41 


36 


31 ± (3) 


50 


0.4 




0.80 < v NN < 0.95 


2.7 ± (0.9) 


1 


29 


2.7 ± (0.9) 
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SS+ 


0.95 < u NN < 0.995 


1.2 ± (0.6) 





34 


1.2 ± (0.6) 


1 


66 




0.995 < u NN < 1.0 


0.6 ± (0.4) 





55 


0.0 ± (0.0) 









0.80 < u NN < 0.95 


8.7 ± (1.6) 


9 


49 


5.7 ± (1.6) 
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SS- 


0.95 < v NN < 0.995 


3.0 ± (1.0) 


4 


36 


3.6 ± (1.0) 


2 


34 




0.995 < v NN < 1.0 


0.9 ± (0.5) 





43 


0.3 ± (0.3) 





70 




0.80 < v NN < 0.95 


169 ± (7) 


169 


50 


73 ± (5) 


64 


19 


FM+ 


0.95 < v NN < 0.995 


55 ± (4) 


43 


9 


19 ± (2) 


18 


49 




0.995 < v NN < 1.0 


20 ± (2) 


20 


48 


3.6 ± (1.0) 


3 


53 



Table 1: The values given in the parentheses are the uncertainties on the mean of 
the background prediction. The Poisson probability for making an observation at 
least as large (or fewer than observed when observed is less than predicted) given the 
predicted background is also shown in the table. The best (blue) and worst (red) 
cases are also shown. 



2.3.3 Results 

The results from the CDF analysis on 3.7 fb _1 are shown in Table [2] and Figure [5j 
6.1±1.8 events were expected in the signal region and 7 events were observed, resulting 
in the observed limits, using 3 NN bins and 5 mass bins, shown in Table |2J 

2.4 DO Analysis with 6.1 fb 1 of Data 

DO uses a similar method as described above. For signal discrimination, variables 
are combined into a Bayesian Neural Network that also uses signal MC and data 
mass sidebands for optimization (Figure E]). Unlike CDF, DO does not have the mass 
resolution to distinguish between and B° s . 
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Channel 


Expected Observed 


B u s Central 


4.0±1.0 3 


B° s Extended 


2.1±0.8 4 


B° d Central 


5.3±1.0 5 


B° d Extended 


2.8±0.8 3 





90% CL 95% CL 




3.6xl(r 8 4.3xl(r s 




6.0xl0- 9 7.6xl0" 9 



Table 2: B® and B d expected and observed number of candidates in central and 
extended mass region (top) as well as resulting observed limits (bottom). 



The DO results are shown in Figure [7J 51±4 background events are expected in 
the highest sensitivity region and 55 events are observed. The expected B° s limit is 
3.8 x 10" 8 at 95% CL observed a limit of 5.1 xlO" 8 . 

2.5 CDF Improvements 

CDF is working on an improved analysis with significantly more (3.2 fb" 1 ) data. 
The forward muon acceptance (left of Figure El) has been increased and the neural 
network has been improved to achieve greater background rejection. The peaking 
background prediction has also been improved. The central-forward (CMU-CMX) 
channel of the analysis has increased in statistics by ~15% resulting in a total increase 
for both channels of ~7%. The neural network has been improved by increasing 
the background rejections (Figure |9]). Additionally the neural network has been 
extensively tested for mass bias (Figure ITU]) . 

The following variables, shown in order of significance to NN separation power, 
are used to discriminate the B® — > fi + fi~ signal from the background: 

Aa: the 3D opening angle between the B® flight direction, p(B), and the direction 
of the decay vertex - estimated as the vector originating at the primary vertex 
and terminating at the muon-pair vertex; 

Iso: the isolation of the candidate B° s defined as, I so = Pt° /(Pt" + YyiPr(^R < 1-0)), 
where the sum is over all tracks within an t] — <p cone of radius R = 1.0, centered 
on p(B); 

|d (/ii)|: impact parameter of muon which has the larger value of the muon pair; 
|d (Bg)|: impact parameter of B® reconstructed from muon pair; 
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Figure 5: Di-muon mass vs NN output distribution for central-central (CMU-CMU) 
and central-forward (CMU-CMX) muons for CDF analysis. 

L2d/o"L 2d : significance of L 2d ; 
Xvtx : X 2 °f t ne vertex; 

L 3D : 3-dimensional vertexes displacement obtained by the vertex fitting; 

Pt(/^2): the transverse momentum of lower momentum daughter muon; 

|do(/i2)|/o"d (/i 2 ) : Significance of smaller |doQf)| of the muon with smaller impact pa- 
rameter, where <Tdo(ji) is estimated uncertainty of \do([i)\] 

\/u\: proper decay length significance (A divided by the uncertainty on A); 

A: the proper decay length; 

|do(/ii)|/o"d ( Ml ): significance of larger |do(^)| of the muon with larger impact param- 
eter, where cr do ^ is estimated uncertainty of \d (fi)\; 

Aa xy : the 2D opening angle between the B° s flight direction, p(B), and the direction 
of the decay vertex - estimated as the vector originating at the primary vertex 
and terminating at the muon-pair vertex; 

|do(/i2)h impact parameter of muon which has the smaller value of the muon pair. 
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Figure 6: Distribution of Bayesian neural network output for signal MC and data 
mass sideband for DO analysis. 



The final NN network consists of 14 input variables. Extensive cross checks were 
done to demonstrate the 14-variable network does not sculpt the di-muon invariant 
mass distribution. Figure ITUl shows the correlation across the M^+^-mass range. No 
significant correlations are seen. 



2.6 Expected Limits for CDF 

A comparison of the previous expected and observed limits is shown in Table [3j The 
expected limit after including the new data set, the increased acceptance, and the 
NN improvements is about B(B° S — > 2 x 10 -8 . 





Expected Observed 


2fb- x 


4.9xl0~ 8 5.8xl0~ 8 




3.3xl0" 8 4.3xl0" 8 


6.9^ 


~ 2 x 10~ 8 



Table 3: CDF expected and observed limits. 
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Figure 7: DO unblinded di-muon mass distributions. The dots with error are the data, 
the solid black line is the expected number of background events, the red dotted- 
dashed lines is the SM signal xlOO, the dashed line is B(D) — > fi + uX and B(D) — > 
fi^vX', and the blue dotted line is B — > n + vD and D — > fi~VX 



3 Conclusion 

FCNC's are powerful probe for New Physics and CDF and DO continue to lead the 
searches in the B sector. CDF is updating their B° sd — > analysis which will 

significantly improve the expected limit for B(B° sd — > fi + fi~). CDF also has the first 
observation of B® — > fi + fi~(j) as well as a A FB measurement that competes and agrees 
with the B-factory results. Shortly after the FPCP conference CDF updated both 
the b — > sji + jjT and B° s d — > /i + /! _ analysis [15j [2j US] . 
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Figure 8: Left: CDF Muon detector layout. The blue region has been added for 
the updated analysis. Right: The fi + fi~K + invariant mass distribution for events 
satisfying the baseline and vertex requirements for the B + — y J/ipK + sample. 
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Figure 10: The top plot shows the NN output distributions for background (black 
histogram) and background trained as "signal" (red histogram). The bottom plot 
shows NN output as a function of di-muon invariant mass. 
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